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We show that the finite difference Backlund formula for the Schrodinger 
Hamiltonians is a particular element of the transformation group on the set 
of Riccati equations considered by two of us in a previous paper. Then, we 
give a group theoretical explanation to the problem of Hamiltonians related 
by a first order differential operator. A generalization of the finite difference 
algorithm relating eigenfunctions of three different Hamiltonians is found, and 
some illustrative examples of the theory are analyzed, finding new potentials 
for which one eigenfunction and its corresponding eigenvalue is exactly known. 
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1. INTRODUCTION 

The factorization of Schrodinger Hamiltonians in terms of differential 
operators of first order plays an important role in the search of quantum 
systems for which the energy spectrum is completely known [1, 2, 3]. It 
is closely related with the existence of intertwining operators [4, 5], Su- 
persymmetric Quantum Mechanics [6] and Darboux transformations [7], 
among other questions. Indeed, most of the exactly solvable potentials 
can be obtained by making use of an appropriate intertwining operator 
transformation. However, nowadays the point is not just to factorize some 
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well-known Schrodinger Hamiltonians [1] but to generate new solvable ones. 
The key to this technique is to intertwine two different Hamiltonians by 
a differential operator A, usually of first order [3, 4, 5, 8, 9, 10] or higher 
[11, 12, 13, 14, 15, 16, 17, 18], provided one of them is solvable. Concern- 
ing iterations of the first order intertwining technique, one of the authors 
(DJF), Hussin and Mielnik [4] have recently used a finite difference algo- 
rithm which provides algebraically the solution of the key Riccati equation 
at a given iteration step in terms of two solutions of the corresponding Ric- 
cati equation at the previous step associated to two different 'factorization 
energies'. This procedure has been successfully applied in order to obtain 
new exactly solvable Hamiltonians departing from the harmonic oscillator 
and Coulomb potentials [4, 5, 19]. 

On the other hand, a number of authors have studied the transformation 
groups related with Riccati equations [20, 21, 22, 23, 24, 25, 26, 27]. In 
the last of these references, it has been used a transformation group on the 
set of Riccati equations in order to analyze their integrability. We think 
that this group theoretic approach could shed light on the abovementioned 
problem of intertwined or related Hamiltonians. A natural question is 
whether there is a relation of this group action with the finite difference 
algorithm. 

We will try to answer these questions in this article, which is organized 
as follows. In Section 2 we recall briefly the problem of two Hamiltonians 
intertwined by a first order differential operator, a technique which is also 
known as A-related Hamiltonians. In Section 3 we relate the Schrodinger 
equations arising in the previous problem with certain Riccati equations 
by making use of the classical Lie theory of reduction by infinitesimal sym- 
metries of differential equations, as applied to homogeneous linear second- 
order equations. In Section 4 we recall the action of a group on the set 
of Riccati equations described in [27]. Using this technique we will prove 
the finite difference theorem [4] in a way alternative to [28]. In Section 5 
we will determine the elements of the transformation group preserving the 
subset of Riccati equations arising from the set of Schrodinger equations 
after applying the reduction process outlined in Section 3. We will find a 
new transformation relating three different Schrodinger equations, which 
represents a generalization of the finite difference Backlund algorithm. In 
Section 6 we find that the problem of A-related Hamiltonians can be ex- 
plained exactly in terms of the transformation group on the set of Riccati 
equations and the reduction procedure of Section 3. Section 7 illustrates the 
use of the new Theorems of Section 5 in the search of potentials for which 
one eigenstate and the corresponding eigenvalue will be exactly known. In 
particular. Examples 7.1, 7.3 and 7.4 will provide potentials essentially dif- 
ferent from the original ones. Finally, we give in Section 8 some conclusions 
and an outlook for future work. 
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2. HAMILTONIANS RELATED BY FIRST-ORDER 
DIFFERENTIAL OPERATORS. 

The simplest way of generating an exactly solvable Hamiltonian H from 
a known one H is just to consider an invertible bounded operator B. with 
bounded inverse, and defining H = BHB~^. This transformed Hamilto- 
nian H has the same spectrum as the starting one H. As^ a generalization 
(see e.g. [3]), we will say that two Hamiltonians H and H are intertwined 
or A-related when AH = HA, where A may have no inverse. In this case, 
if V' is an eigenvector of H corresponding to the eigenvalue E and A^f) ^ 0, 
at least formally Atp is also an eigenvector of H corresponding to the same 
eigenvalue E. 

If ^ is a first order differential operator, 

A=-^ + W{x), and A^=-± + W{x), (1) 
then the relation AH = HA, with 

F = -^ + y(.), H = -^ + y(.), (2) 

lead to 

V = -2W' + V, W{V -V) = -W" - V . 

Taking into account the first equation, the second becomes 2WW' = W" + 
V' , which can easily be integrated giving 

V==W'^-W' + e, (3) 

and then, 

V = W'^ + W' + €, (4) 

where e is an integration constant. The important point here is that H 
and H, given by (2), are related by a first order differential operator A, 
given by (1), if and only if there exist a constant e and a function W such 
that the pair of Riccati equations (3) and (4) are satisfied simultaneously. 
Moreover, this means that both Hamiltonians can be factorized as 

H = + e , H = AA'i + e. (5) 

Adding and subtracting equations (3) and (4) we obtain the equivalent 
pair which relates V and V 

V-e = -{V~e) + 2W'^, (6) 

V = V + 2W'. (7) 
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The function W satisfying these equations is usually called superpotential, 
the constant e Js the factorization energy or factorization constant and V 
and V (resp. H and H) are said to be partner potentials (resp. Hamilto- 
nians). 

Notice that the initial solvable Hamiltonian can indistinctly be chosen 
as H or H. In both cases the point will be to find a solution W of the 
corresponding Riccati equation (3) or (4) for a specific factorization energy 
e. Prom this solution the expression for the (possibly) new potential follows 
immediately from (7). 



3. DILATION SYMMETRY AND REDUCTION OF A 
LINEAR SECOND-ORDER DIFFERENTIAL EQUATION 

Let us briefly recall a well-known method of relating a homogeneous 
linear second-order differential equation to a Riccati equation, which can 
be regarded as an application of the classical Lie theory of infinitesimal 
symmetries of diff'erential equations. Its importance will become clear when 

applying the method to timc-indcpcndcnt Schrodingcr equations. 
The homogeneous linear second-order differential equation 

^+K-)-+c(.). = 0, (8) 

admits as an infinitesimal symmetry the vector field X = zd/dz generat- 
ing dilations (sec e.g. [29]) in the variable z, which is defined for z 7^ 0. 
According to the Lie theory of infinitesimal symmetries of differential equa- 
tions, we should change the coordinate zto a. new one, u = fiz), such that 
the vector field X = zd/dz becomes a translation generator X = d/du 
in the new variable. This change is determined by the equation Xu = 1, 
which leads to u = log \z\, i.e. \z\ = e". In both cases of regions with z > 
or z <0 we have 

dz du , d'^z f d.u\'^ d'^u 

= z— , and -tt; — z\ — \ +z- 



dx dx ' dx"^ \dx J dx"^ ' 

so the equation (8) becomes 

d'^u , , .du f d.u\'^ , . 

d^^'^'^Tx^[Tx) = 

As the unknown function u does not appear in the preceding equation but 
just its derivative, we can lower the order by introducing the new variable 
w = du/dx. We arrive to the following Riccati equation for w 

—— = —w^ — b(x)w — c(x) . (9) 
dx 
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Notice that from dz/dx = zdu/dx and the definition of w we have 

1 dz 

w^- — . (10) 

z dx 

The second order differential equation (8) is equivalent to the set of (9) 
and (10), because given a function w satisfying (9), the function z defined 

(up to a factor) by (10), i.e. z(.t) = cxp 

(/''w(C)dC), satisfies (8). Wc 
could have followed a similar pattern straightening out the vector field in 
the opposite sense, that is, by imposing Xu = —1. This would have lead 
to zt = — log 1 2: 1, or \z\ = e~". Now, in either case of 2; > or ^ < we have 
dz/dx = —zdu/dx and d'^z/dx^ = z{du/dx)'^ — zd'^u/dx^, so we finally 
obtain the Riccati equation 

^ = - b{x)w + c{x) , (11) 



where now 



du 1 dz 
dx z dx 



We will distinguish in what follows between these two alternatives of 
reduction of (8) by means of a subscript + or — in the corresponding 
functions, respectively. Wc remark that both are defined locally, that is, in 
open intervals where z has a constant sign. 

Let us apply these ideas to the particular case of the one- dimensional 
time-independent Schrodinger equation 

-^ + {V{x)-e)cl> = 0, (13) 

where V{x) is the potential and e is some specific energy eigenvalue. As 
explained before, we can reduce (13) either to the pair 

W[ = -Wl + iV{x)-e), WV = ^£, (14) 

or, alternatively, to the pair 

WL=Wl-(V(x)-e), W- = -\^. (15) 

dx 



The Riccati equations appearing in these pairs resemble those appearing in 
Section 1, namely equations (3) and (4), but in the systems (14) and (15) 



6 



CARINENA, FERNANDEZ AND RAMOS 



the unknown functions W+ and W- are related by VF+ = — W-, while in 
both (3) and (4) the unknown W is the same function. 

However, the previous remark will be useful in the interpretation of equa- 
tions (3) and (4). We can rewrite them as 

W' = W^-{V{x)-e), (16) 
W' = -W^ + {V{x) - e) . (17) 

Then, we can regard equation (16) (rcsp. equation (17)) as coming from a 
Schrodinger-type equation like (13) (resp. like — (P4>/dx'^ + {V{x) — e)(j) = 0) 
by means of, respectively, the changes 

^ (j) dx ^ ^ ^ dx ' ^^^^ 

so the two "eigenfunctions" (p and (p of the mentioned Schrodinger-type 
equations are related by (pep = Const. Of course, the changes (18) are 
defined locally, i.e. in common open intervals of the domains of (f> and (f> 
determined by two consecutive zeros of </> or 0, or maybe by a zero and a 
boundary of the domain of the problem. Note that there is no reason why 
they should provide functions (p, (j) defined in the same way in the entire 
domain of W , but in general they will be defined interval-wise. Moreover, 
if we choose the function W of the two operators A and defined in (1) 
as given by (18) it holds Acp = Q and cj) = 0. 

We have seen that the Riccati equations (16) and (17) correspond by 
means of the changes (18) to two Schrodinger-type equations which in turn 
are equivalent to 

H<i) = e(i), H^=e^, (19) 

where H and H are given by (2). Then, it is equivalent to say that H 
and H are ^-related, with associated constant e, to say that the pair of 
functions (p and (p, which satisfy = Const., are respective eigenfunctions 
with eigenvalue e of the Hamiltonians H and H. Each of these facts imply 
that both Hamiltonians can be factorized as in (5). We finally remark 
again that these factorizations make sense only locally, i.e. in common 
open intervals where (j) and (p are defined. ^ 
A special case where all becomes globally defined arises when (p ov (p 
is the ground state wave function of its respective Hamiltonian, having 
then no zeros in the entire domain of the problem. Then, remembering 
the equation (7) it is clear the relation of what we have just exposed with 
the Darboux transformations in the context of one-dimensional or Super- 
symmetric Quantum Mechanics, as it is shown, for example, in [7, pp. 7, 
24]. 
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Moreover, for e below the ground state energy of H (resp. H) it is 
sometimes possible to find a non-normalizable eigenfunction (resp. (p) 
of H (resp. H) without zeros, leading to physically interesting potentials 
[4, 5]. However, as this procedure is slightly involved, we shall restrict for 
simplicity the discussions along this jirticle to the previous case of e equal 
to the ground state energy of H or H. 

4. TRANSFORMATION GROUP ON THE SET OF RICCATI 

EQUATIONS 

In this section we introduce the technique of relating different Riccati 
equations by means of a group action. The procedure, of key importance 
in what follows, was described in [27] in order to give a group theoretical 
explanation of some properties of the Riccati equation. 

We recall that a general Riccati equation is the non-linear first order 
differential equation 

= 02 {x) y2 (x) + ai {x) y{x) + ao {x) , (20) 

ax 

where the coefficients ai{x), i — 0, 1, 2, are smooth functions of the inde- 
pendent variable x. A particular Riccati equation, determined by a specific 
election of these coefficient functions, can be regarded as a curve in M'', i.e., 
an element of Map(M, M^). 

We can transform every function y{x) in the extended real line M by an 
element of the group of smooth SL(2, M) valued curves Map(M, SL{2, M)), 
which from now on will be denoted by means of the action O : ^ x 
Map(]R, 1) ^ Map(IR, 1) defined as follows [27]: 

9(A(.),,(.))= ";-f;-'-;+f(-| , if.(.)^-M, (21) 

e(A(.),oo) = ^, e(^(x),-M)=oo, (22) 

when 

It is easy to check that the Riccati equation (20) transforms under the 

generic transformation y{x) = 0{A{x). y{x)) into a new Riccati equation 
for y{x) with coefficients ai{x), i = 0, 1, 2. The relation amongst the new 
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\ao) 






and the old coefficients is given by 

02 = S'^ a2 — ai + 7^ ao + ^5' — 5j' , (24) 
ai = -2i35a2 + {aS + ai - 2 07 oq + 5a' - a5' + fi^' - 7/?' , (25) 
So = 02 — a/3 fll + a^ flo + a/?' — /3a' , (26) 

where the prime denotes derivative respect to x. We would like to remark 
that similar transformations, using GL{2, R) instead of SL{2, M) have been 
used in [30, 31], also in connection with transformations of the Riccati 
equation. However, it seems that none of them noticed that these trans- 
formations constitute indeed a group action, as we will see next. 
The previous relation can be written in a matrix form 

52 -6'y 

2/3(5 aS + (3"f - 
/32 -a/3 

7(5' - sy 

+ 1 (5a' - a6' + fiy - 7/3' | . (27) 
a/3' - /3a' 

If we define the maps 

(52 -,57 72 
B{A) = ( -2/35 a(5 + /37 -2a7 | , (28) 
/32 -a/3 a2 

7^' - 6'y' 

e{A) = I 6a' -aS' + I3Y -^13' | , (29) 
a/3' - /3a' 

where ^ G ^, it can be easily checked that B is a linear representation 
of the group Q on Map(]R, M.^). The restriction of B to the subgroup 
of constant SL{2, ]R)-valucd curves is nothing but the adjoint action of 
SL{2, M) on its Lie algebra [27]. On the other hand, it can be proved that 
9 defines an 1-cocycle for B, which means that (see [27, 32]) 

9{A,A2) = B{A^){e{A2)) + 9{A^) , \/Ai,A2Gg, (30) 

and moreover 9 is not an 1-coboundary for B. Consequently, see e.g [32], 
the expression (27) defines an affinc action of Q on the set of general Riccati 
equations, which in turn can be identified with the set of curves on the Lie 
algebra sl(2, R). 

In practical terms, all this means that the composition of two transfor- 
mations of type (27) on the set of Riccati equations is exactly the same as 
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making only one transformation in the following way. If Ta denotes the 
transformation of type (27) associated with A G G, then it holds 

Ta, o Ta, = Ta,a, , W Ai, A2 G G . (31) 

It is possible to generalize this afBne action in more general situations, 
when an arbitrary finite dimensional Lie group is involved [33] . 

In short, if y{x) is a solution of the equation (20) and we transform it by 
means ofy{x) = Q{A{x),y{x)), then y{x) will be a solution of the Riccati 
equation with coefficient functions (24), (25) and (26). By means of this 
technique we will prove the following result (the notation has been adapted 
to the one used along this article). 

Theorem 4.1 (Finite difference Backlund algorithm [4, 28, 34, 35]). 

Let Wk{x), u!i{x) be two solutions of the Riccati equations w' + w'^ = 
V{x) — Ck and w' + w"^ = V{x) — ei, respectively, where < e;. Then the 
function Wki {x) defined by 

Wki{x) = -Wk{x) -p\ — ^Vt' (^2) 

Wk{x) - wi{x) 

is a solution of the Riccati equation w' + = V{x) — 2 w^(x) — e;. 



Proof. The function 'Wi{x) satisfies the Riccati equation w' + w"^ = 
V{x) — e; by hypothesis. We transform it by means of the element Ao{x) 
of G given by 

= + ,33, 

where h{x) and a are a function with the same domain as wi{x) and a 
positive constant, respectively. Notice that Af){x) G G since its determinant 
is always one, for all x in the domain of h{x). According to (21) and (22) 
we compute 

. , ^ . h(x)'Wi(x) — h'^(x) + a , " 

e{Ao{x),wi{x)) = ^ '\\' ^ = -h{x) + 



h{x) - wi{x) h{x) - wi{x) ' 

This is a solution of the Riccati equation with coefficient functions given 
by (24), (25) and (26), with matrix elements a{x) = 5{x) = h{x)/^, 
/3{x) = (— /i^(x) + a)/y/a, j{x) = —1/y/a and coefficients of the initial 
Riccati equation 02(2;) = —1, ai{x) = and ao{x) = V{x) — e;. Simply 
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performing the operations, we find 

Mx) = ^{-h^{x)-h'ix) + V{x)-ei+a}-l, 

ai{x) = ^^^{-h''{x)-h'{x)+V{x)-€i+a}, 

ao{x) = ^^{-h\x)-h'{x)+V{x)-ei+a} 
+h^{x) + h'{x) - 2h'{x) - a . 

Therefore, if the function h{x) satisfies the Riccati equation + w' = 
V{x) — €k, with cfe = e; — a, and we rename it as h{x) = Wk{x), the new 
coefficients reduce to a2{x) = —1, ai{x) = and ao{x) = V{x) — 2 w'f.{x) — 
ei. I 

Let us note that in [4] the proof of the Theorem 4.1 was just sketched. 
In addition, there exists an alternative proof; see, e.g., Mielnik, Nieto and 
Rosas-Ortiz [28]. 

5. GROUP ELEMENTS PRESERVING THE SUBSET OF 
RICCATI EQUATIONS w' + w"^ =^ V{x) - e 

We have just seen how the transformation group defined on the set of 
Riccati equations provides a direct proof of Theorem 4.1. It relates one 

solution wi{x) of the initial Riccati equation w' + w'^ = V{x) — e; with one 
solution Wki{x) of the final Riccati equation w' + uP' = V{x) — 2ui'f.{x) — 
€[ by using a solution Wkix) of the intermediate Riccati equation w' + 
w'^ = V{x) — Efc. These three Riccati equations can be obtained from 
another three Schrodinger-like equations by means of one of the reduction 
possibilities explained in Section 3. Moreover, those associated with the 
initial and intermediate Riccati equations, namely —ijj" + {V{x) — €i)ijj = 
and —ip" + {V{x) — ek)^ = 0, can be seen as the arising eigenvalue 
equations for the two energies e;, tk of the same potential V{x), meanwhile 
the final Riccati equation can be assocdatcxi to the eigenvalue equation for 
the potential V{x) — 2w'f.{x) with eigenvalue e;. 

Then, we are naturally led to the question of which are the most general 
elements of Q preserving the subset of Riccati equations characterized by 
having the coefficients 02(2;) = —1, aiix) = and a^ix) equal to some 
function, which we will write as an expression of the form V{x) — e. 

The Riccati equation we will start from is 



w' = -w"^ + V{x) - e. 



(34) 
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which according to (20) has the coefficients a2{x) = —1, ai{x) = and 
ao(a;) = V{x) — e. The condition for obtaining a final Riccati equation in 
the mentioned subset is 

-2j35 a5 + (3-f -2a-f 
7<5' - .57' 

+ 1 Sa'- aS' + - 7/3' 1 ' (^5) 
a/3' - I3a' 

for an A{x) G Q given by (23) to be determined, and where V{x) — e will 
be in general different to V{x) — e. Therefore, the elements of the subset of 
G we are trying to characterize will not necessarily form a subgroup. The 
matrix equation (35) is equivalent to three scalar equations 

-1 = -S^ +"f^ {V{x) - e) +"f6' ~ 6"f' , (36) 
= 2/3(5-2 a7 {V{x) - e) + 6a' - aS' + Pi - 7/3' , (37) 
V{x) - e = -p'^+a^ (y(x) - e) + a/3' - pa' . (38) 

Differentiating det A(a:) = a{x)5{x) — P{x)^{x) = 1 we have as well 

a' 5 + 5' a - iP -p'^ = Q. (39) 

Out of these four equations, (36), (37) and (39) will give conditions on the 
matrix elements a, /3, 7, 5 and their derivatives such that the preserving 
condition be satisfied. The remaining (38) will define V{x) — e in terms of 
all the other functions, including V{x) — e. 

After taking the sum and the difference of (37) and (39) it follows 

{V{x) - e)a^ = — + a/3' , (40) 

7 7 

{V{x)-eh^ = ^-f + ^-lS'. (41) 
Substituting them into (36) and (38) gives 

-1 = -6^ + ^ + ^-5i, 
a a 

TT/ X ^9 ctPS 5aa' „ , 

V{x) - e = -P^ + — + pa' . 

7 7 
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Multiplying the first of these equations by a and the second by 7, and 
using the fact that — /J7 = 1, we arrive to 

a = J + 7', (42) 
(F(a;) - e)7 = /3 + a' . (43) 

Substituting (42) into (41) yields 

{V{x) - e)j = I3-S'. (44) 

We have two relations amongst the functions a, /?, 7 and 5, namely (42) 
and the determinant condition, so we can express these matrix elements in 
terms of only two of them and their derivatives. Then we have a = i5 + 7' 
and P = {5{5 + 7') — l)/7. Using moreover the fact that 

7 V7y 7^ 

the equation (44) becomes 

5\' ( SV . 1 



7/ V 7y 7 
so the new function v defined as u = — 5/7 must satisfy the Riccati equation 

v' + = V{x) + ^-e. (45) 

Now, substituting in (43) the expressions of f3 and a' in terms of d, 7 and 
their derivatives, using the definition of v and the equation (45) gives 

V{x) - e = V{x) - 2 u + v'^ + ^ - e . 

It only remains to find the expression of the function solution of the final 
Riccati ciquation, in terms of w and v. The 5^(2, R) -valued curve used for 
the transformation can be written as 

Co{x)=^ "^-"^y^-^^ , (46) 

so the desired function is 

—vw + wj' /j — l/j"^ +v'^ — vj' /j 



w{x) = @{Co{x),w{x)) 



1 - - y 

= -t; - + -L . (47) 

w — V 7 
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In short, we have just proved the foUowing Theorem: 
Theorem 5.1. Let w{x) be a solution of the Riccati equation 

w' + = V{x) - e (48) 

for some function V{x) and some constant e, and 7(2;) a never vanishing 
differentiable function defined on the domain ofV{x). Ifv{x) is a solution 
of the Riccati equation 

v' + = V{x) + ^^-e, (49) 
7^ (a;) 

such that is defined in the same domain as w{x) and w{x) — v{x) does not 
vanish, then the function w{x) defined by 

Mx) = -vix) - + (50) 

w{x) — v{x) j{x) 

is a solution of the Riccati equation 

W' = V{x)-2(j-v + v'^ +^-e. (51) 



Needless to say, the coefficients of the final equation can be calculated 
directly by using (24), (25), (26) and taking into account (46), (48) and 
(49). 

Corollary 5.1. The Theorem 4-1 is a particular case of Theorem 5.1. 

Proof. It is sufficient to choose in Theorem 5.1 w{x) = wi{x), v{x) = 
Wk{x), € = ei and 7 = I/Vq - efe, with Ck < ei. | 

Theorem 5.1 has a counterpart for linear second-order differential equa- 
tions of Schrodinger type, which will be in turn of direct interest in physical 
applications. The key is to use in a inverse way the reduction procedure 
outlined in Section 3. 

Consider the solution w of the Riccati equation (48). We can define 
(locally and up to a non-vanishing multiplicative constant) the new variable 
4'w as 

Mx)^exp( [ wiOd^) , (52) 
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which will satisfy 

for the specific constant e. Analogously, by considering a solution v of the 
Riccati equation (49) we can define (locally etc.) (j)y as 

</),(ar)=exp(^y%(Ode), (53) 

which will satisfy 

-K+(nx) + ;^-e)</..=0, 

for the same specific constant e. Then the function w defined by (50) will 
satisfy the Riccati equation (51). We could define as well (locally etc.) the 
new function as 



(54) 



(f>w{x) =exp(^J w{^) , 
which in turn will satisfy 

What has to be done now is to relate the function (j)yj with (p^ and (py, 
taking into account the relation amongst W, w and v. 

Proposition 5.1. Let w, v, W be the functions for which the Theo- 
rem 5.1 holds, and (j)w, (pv, <t>w the ones defined by (52), (53) and (54), 
respectively. Then we have 

^=^, ^=,, J=^, (55) 

(Pw <Pv 9w 



and it holds 



dx 

up to a non-vanishing multiplicative constant. 



<fe = 7{-A + ^U^, (56) 
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Proof. The first assertion is imniediate. As a consequence, we have 
^(~lx + ^)'^«' = '^(^ ~ Taking the logarithmic derivative 

{jiv - w)(p^y ^ y + ^' ~ ^' _^ ^ 

7(1; - 7 v-w 4>u, 

= — H h — ^-^ hw = — -tw-'yH — \-w 

7 v~w v—w 7 V— w 



7 u - w ^ ' 

where it has been used the equations (48), (49) and (50). | 
With the previous results we have the following: 

Theorem 5.2. Let (t>w{x) be a solution of the homogeneous linear second 
order differential equation 

-<j>'i + {V{x)-e)<i>^ = 0, (57) 

for some specific function V{x) and constant e, and 7(0;) a never vanish- 
ing differentiable function defined on the domain of V{x) . If the function 
(j)y {x) ^ (pw {x) is a solution of the equation 



-4,'^+\^V{x) + ^-ej^y=0, (58) 

defined in the same dom.a,in as (p^ix), then the function (j>w{x) defined (up 
to a non-vanishing multiplicative constant) by 

'fe = 7f-:^ + T'V^, (59) 



dx 



satisfies the new equation 



V{x)~2{^v + v'] +^-e}(j>^^0, (60) 



where the function v{x) is defined (locally) as (p'y/(l>v = v. 



Note that Theorems 5.1 and 5.2 are invariant under the change of sign 
of 7. 
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6. FINITE DIFFERENCE ALGORITHM AND 
INTERTWINED HAMILTONIANS FROM A GROUP 
THEORETICAL VIEWPOINT 

We have already said that the finite difference algorithm, based on the 
Theorem 4.1, appeared in [4] when the authors wanted to iterate the stan- 
dard first order intertwining technique. This idea has been kept also in 
subseqiient works [5, 19], and in all of these articles the algorithm has 
been shown to be of use for obtaining new exactly solvable Hamiltonians. 
Moreover, the proof of Theorem 4.1 given recently by Mielnik, Nieto and 
Rosas-Ortiz, alternative to that which has been given here, still relies on the 
idea of iteration of the intertwining technique, see [28, Sec. 2] for details. 

On the other hand, we have given a direct proof of Theorem 4.1 by 
making use of the action of Q on the set of Riccati equations, and we 
wonder whether it is possible to establish a further relation between this 
transformation group and the (maybe iterated) intertwining technique. 

The important result, which we show next, is the following. By using 
properly the finite difference algorithm just once, jointly with the reduc- 
tion procedure described in Section 3, it is possible to explain from a group 
theoretical viewpoint the usual problem of ^-related or intertwined Hamil- 
tonians. 

With this aim, let us consider two Hamiltonians 

= + ^^ = -^ + ^^^'^)' 

which by hypothesis are Ai-related, i.e. AiHi = HqAi and HiA\ = AIHq, 
where 

Ar = ^+-r, 4=-^+-i, (62) 

and wi is a function to be determined. 

Assume that Hq is an exactly solvable Hamiltonian for which we know 
a complete set of square-integrable eigenfunctions tpn^ with respective en- 
ergies En, n = 0, 1, 2, We have seen in Section 3 that, in particular, 

Vo{x)-Eo = wl{x,Eo)+w[{x,Eo), (63) 
Vi{x)-Eo = wl{x,Eo)-w[{x,Eo), (64) 

or equivalently 



Vo{x)-Eo = -iVAx)-Eo) + 2wl{x,Eo), (65) 
Vi{x) = Vo{x)-2w[{x,Eo). (66) 
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where we have chosen wi{x, Eq) as 

t«i(a;,Eo)= V-r/^r- (67) 

Up to a non- vanishing multiplicative constant, we define the function tp^'^ 
as ■^0^' = We have as well 

wi{x,Eo) = -4'^'/4'K (68) 
Then, both Hamiltonians factorize as 

Ho = {Eo)A\ (Eo) +Eo, H^= A\ {Eq)A^ (Eo) + E^ . (69) 

We have made explicit Eq in the function Wi and, as a consequence, in the 
operators Ai and . However, it should be considered as a label reminding 
the factorization we are working with rather than as a functional depen- 
dence. From (67) and (68) we have A\{Eq)iI)^^^ = and Ai{Eq)^''^'^ = 0; 
as a result Hiipj^^ = Eoipj^^ and Hqi/jI^^ = Eoijjj^K As has no zeros 
in the domain of Vq{x), all the functions defined in this Section will be 
globally defined provided that such a domain is connected. 

The equation (66) relates the new potential Vi (x) and the old one Vq {x) . 
As it is well known, due to the Ai(£Jo)~relationship of the Hamiltonians 
Hq and Hi, the normalized eigenfimctions of Hi can be obtained trans- 
forming those of Hq by means of the operator A\{Eo) except Vo^'*' since 
A\(Eo)'tj}Q^^ = 0. In fact, a simple calculation shows that the functions 

satisfy 

ifi^i) = and (V;i^),^W)=<5„™, (71) 

for all n, m = 1, 2, 3, . . . , provided that the functions ipn'^ are normalized. 

Although the function ■^q^'' satisfies Hiipo'^ = £'oV'o^\ it does not corre- 
spond to a physical state of Hi since it is not normalizable, which means 
that Eq will not belong to the spectrum of Hi . For this reason the Hamil- 
tonians Hi and Hq are said to be quasi-isospectral. 

Let us formulate now these results in terms of the transformation group 
on the set of Riccati equations introduced in Section 4. By hypothesis we 
have 

ifoV'i°' = n = 0, 1, 2, .... (72) 
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As Hq is given by (61), the spectral equation (72) can be written as the 
following set of equations 

-V'i°^"+(Fo(x)-i?„)Vf =0, n = 0, 1, 2, ... (73) 

We introduce the new functions 

,(0)/ 

w^{Er,) = ^, n = 0, 1, 2, (74) 

where the dependence on x has been omitted for brevity. As we know 
from Section 3, these transformations will be defined locally, i.e., for each 
n the domain of wi(En) will be the union of the open intervals contained 
between two consecutive zeros of ipn '' or maybe a zero and one boundary 
of the domain of Vq{x). In particular, w-i{Eo) is defined globally in the 
entire domain of Vo(x) as ipQ^"^ have no zeros there. Therefore, the set of 
equations (73) reads in the new variables as the set 

w[{En) + wl{En) = Vo{x) -En, n = 0, 1, 2, . . . , 

that is, the functions wi{En) are respective solutions of the Riccati equa- 
tions 

w' + up = Va{x) - E„ , n = 0, 1. 2, .... (75) 

Let us apply the Theorem 4.1 to this situation. We act on the set 
of all equations (75) but one by means of suitable group elements of Q. 
These 5L(2, R)-valued curves are constructed by means of the solution of 
the equation of the set (75) which is to be set aside. In order to avoid 
singularities, this solution should be the one with n = 0. 

The mentioned elements of Q are analogous to the one used in the proof 
of Theorem 4.1. They turn out to be 

V-E„ - Eo V -1 ^i(^o) ; ' 

(76) 

We define the new functions wi{En) by 

- iir \ r^iTi rir w Wi{Eo)wi{En) - wI{Eq) + En - Eq 
Wl{En) = e{Bn,Wl{En)) = ^ '- — 

Wi(Eq) - Wi{En) 

= -Wi{Eo) i^\^^"ir, . , n=l, 2, .... (77) 

By Theorem 4.1 these functions satisfy, respectively, the new Riccati equa- 
tions 

w' +w'^ = Vi{x)- En, n = l, 2, (78) 
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where Vi{x) — Vo{x) — 2w[{Eo)- We can define (locally etc.) the new set 
of functions (j)n \ for n = 1, 2, . . . , by 

.^W(x)=exp(^y wi{^,En)d^y n=l, 2, (79) 

which therefore satisfy, respectively, a linear second-order differential equa- 
tion of the set 

-ct>" + (Viix) - Er,)(j, = , n = l, 2, .... (80) 

Then, <pn ^ are eigenfunctions of the Hamiltonian Hi = —-j^ + Viix) with 
associated eigenvalues En, for n = 1, 2, .... As a consequence, they can be 
written as the hnear combinations (j)n'' = ipn^ + A„ ipn^x) , for 

all n = 1, 2, . . . , up to non- vanishing constant factors, where A„ are still 
unknown constants and it has been used the well-known Liouville formula 
for finding the second linearly independent solution of each equation of (80) 
starting from xp'^^"' . 

Now, the important point is that each of the functions 0n ■* turns out to 
be the same as ipn^ , up to a non- vanishing constant factor, i.e. the previous 
constants A„ are all zero. In fact, the logarithmic derivative of ipn^ is 

^ A(4(i,„)^(°)) _ i^Mi"^' + ^4°^) 

4(i^o)4°^ _^(o)' + vr^(o) 

I (0)" , V'o"'" ; (") f V^o"'' ^ / (0) , V-^,"'' / (0)' 

Yn ^ . (0) Yn 
Wo 

Taking common factor ipn^ in both numerator and denominator, using the 
relations 

; (0)" 

^=Vo{x)-En, n = 0, 1, 2, 
and the definitions (74), we arrive to 

-^ = -Wi{Eo) — — ——= wi{En) = , (81) 

for n = 1, 2, . . . , therefore ^j^^'' and (pl^^ must be proportional. It is also 
clear that these equations hold interval- wise. 
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Now, as far as ■^q ' is concerned, it is clear that Theorem 4.1 would 
make no sense for Wk{x) = wi{x) and Ck = ei- In a similar way, we cannot 
put En = E() in (76): the normalizing factors 1 / \/ £"„ — Eq , which were 
introduced in order to get S'L(2, R)-valued curves, would no longer make 
sense because these matrices, after dropping out such factors would have 
zero determinant. This means that one cannot use a transformation of type 
(77) with En = Eq for the function wi{Eo) itself. However, the associated 
function to ipo^^ at the Riccati level is just given by (68), i.e. the new 
function Wi{Eo) = —'w-i{Eo) satisfies an equation of type (78) for n = 0, 
which is precisely (64). In short, the equation A\{Eq)'(P^^ = is translated 
at the Riccati level into the fact that wi{Eo) cannot be transformed in the 
sense mentioned above. Conversely, it is not possible to write wi{Eo) = 
e{Ba,wi{Eo)) for Bo G G. 

We have just explained the problem of two Ai{Eo)-Ye\a.ted Hamiltonians 
by means of Theorem 4.1, which in turn is a particular case of Theorem 5.1. 
For the sake of completeness, let us show briefly how Theorem 5.2 applies 
to the same problem. Consider the set of equations (73), where we retain 
again the one with n = 0, which will play the role of equation (58). All 
the others will play the role of equation (57). For each n = 1, 2, . . . , the 
function 7 must be defined by 

—Eq = —En H — ^ . 

T 

Thus, we can choose 7 = l/\JEn — Eq. According to (59) and (60) the 
functions 

satisfy, respectively, 

+ iVo{x)-2w[{Eo)-En)^ = 0, n=l, 2, ... , 

where it has been used wi(So) = tp'^^' / '^q'^ . In this way we have recovered 
the normalized eigenfunctions (70) associated to the new potential Vi {x) = 

Vo{x)-2w[{x,Eo). 

7. ILLUSTRATIVE EXAMPLES 

In this Section we will apply Theorem 5.2 to some cases where j{x) is not 
a constant, what provides a more general situation than the usual intertwin- 
ing technique. However, note that with this method we will be able to find 
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potentials for which one eigenfunction and its corresponding eigenvalue are 
exactly known. We will use a slight generalization of two well-known types 
of potentials, namely the radial oscillator and Coulomb potentials, which 
consists of taking the most general intervals of the appearing parameters 
such that it is possible to find square- integrable eigenfunctions. 

7.1. Radial oscillator- like potentials 

Let us consider the family of potentials 

Vi,t{x) = — + , (83) 

where x G (0, oo) and I, b are real parameters. Each member can be re- 
garded as being part of a pair of shape invariant partner potentials, with 
associated transformation law I l + l for the parameter I [1, 9]. This fact 
allows to find the eigenvalues and the corresponding eigenfunctions, even 
normalized, in an algebraic way. The key is to find functions in the kernel 
of the first order differential operator d/dx — {l + l)/x + fox/2 which be nor- 
malizable with respect to the norm induced by the standard scalar product 
defined in £^(0, oo). That will provide the groimd state eigenfunction. The 
eigenfunctions of the excited states are obtained from the iterated action 
of —d/dx — {I + l)/x + fox/2, with appropriate I at each step, times some 
suitable factors, on the ground state eigenfunction. However, the point 
is that with this procedure, one obtains the boundary conditions of the 
eigenfunctions as a consequence rather than being a priori requirements. 
The result for this family of potentials are the normalized eigenfunctions 
(up to a modulus one factor) [36] 



where k = 0, 1,2, . . . , for / > —3/2 and fo > 0. The notation L^{u) means 
the Laguerrc polynomial of degree n and parameter a in the variable u. 

Note that in the interval / G (—3/2, —1) these eigenfunctions go to infinity 
as X tends to zero, contrary to the usual requirement of going to zero, in 
spite of the fact of being squarc-intcgrablc. 

The problem of the quantum-mechanical motion of a particle in a poten- 
tial on the half line (0, oo) has been carefully studied in [37]. It has been 
shown there that the operator H = —■^ + V{x), with domain C^(0, oo) of 
differentiable functions with compact support in (O. oo), V being a contin- 
uous real- valued function on (O.oo), is a symmetric operator and that it is 
essentially sclf-adjoint if and only if V{x) is in the limit point case in both 
zero and infinity [37, Theorem X.7]. In the case we have in hand, what 
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happens is that the potentials of the family (83) lead to essentially self- 
adjoint Hamiltonians —(P/dx^ + Vifi{x) for the range 6 > and I > —3/2, 
with different sclf-adjoint extensions in each of the intervals I G (—3/2, —1) 
and I e (— l,oo), the first including eigenfunctions which do not necessar- 
ily go to zero as x — > 0. We will see that one eigenfunction arising when 
/ G (—3/2,-1) provides an interesting application, for the family of po- 
tentials (83), of our new method generalizing the first order intertwining 
technique. 

In both cases, the corresponding eigenvalues to the eigenfunctions (84) 
for the potentials (83) are [36]: 

E^^^ = b(2k + l + ^] , fc = 0, 1, 2, .... (85) 



2, 

If 6 = 2 these eigenvalues reduce to those of [38]. Compare also with the 
eigenfunctions and eigenvalues given in [15, pp. 391-2]. 

Example 7.1. Let us consider the following variant of the family of 
potentials (83) 

where x G (0,oo), I > —3/2 and b > 0. The normalized eigenfunctions 
are given again by (84), with the same peculiarities, but the corresponding 
eigenvalues are now: 

E^^^ = 2bk, fc = 0, 1, 2, .... (87) 

We would like to apply Theorem 5.2, by using two potentials of the family 
(86) for different specific values of I. The difference between them should 
be a positive function in (0, oo) in order to define appropriately 7(2;) as 
required by the Theorem. We have 

/ 21 

Vifiix) - Vi+r,b{x) =r\b 

where r is a positive integer. Since 6 > 0, the condition for the right hand 
side to be always positive is that 2Z + 1 + r < 0. We can find a solution 
if r = 1, since then it should happen 2/ + 2 < or equivalently / < — 1. 
For r = 2, 3, . . . , we would find I < —3/2, which is incompatible with the 
range of I. Then, we have to choose r = l, — 3/2<Z<— 1 and therefore 
an appropriate function 7(0;) is 



lhb{x) 



2l + o\ -1/2 
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We will transform an eigenstate of Vi+i^b{x) by making use of the eigenstate 
o{Vifi{x) with the same energy, i.e. with the same k. Consider the functions 



1 



one of which will be used to find the final potential according to (60) 



Due to the presence of the Laguerre polynomials in (84), Ck''{x) has k 



zeros in (0,oo) and therefore v''il^{x), as well as the final potential, have k 

singularities in the same interval. In order to avoid them, wc choose k = 0. 
Summarizing, we transform the eigenfunction Co^^'^i^) obeying 

-^^^^ + Vt+MC''\x) = , (88) 

by means of the solution (q^{x) of a equation similar to (88) but with / 
instead of ^ + 1. Since I G (—3/2, —1), both of the original eigenfunction 
Co^^' (x) and the intermediate one (x) are square-integrable, but this 
last goes to infinity when x 0. 

After some calculations, the image potential becomes 

T/im/ N / \ rif^l'lji l,b , ^'-'o''\ , "T/.i 

V7/,b dx J -fub 

Px"^ (/ + !)(/ + 2) / 3\ 6b{l + l) 



+ ' ^ -b l + - 



2 J (&a;2 - 2(Z + 1))2 
for which we obtain the eigenstate with zero energy 

Vo (x) = 7i,b{x){ ^ +^^0 Co (x) 



/ 1,1+5/2 ^l+2^-bx^/i 

" V 2'+3/2r(/ + 5/2) ^bx^ - 2{l +"fy ' 

as can be checked by direct calculation. Notice that hx'^ — 2(? + 1) > 
always since / < — 1 and 6 > 0, and therefore VJ*™(x) and rj^ {x) are defined 

in the whole interval (0, oo). Moreover, 'I]q^{x) has no zeros, and it tends 
to zero when x goes to and to oo fast enough to give a square-integrable 
eigenfunction. In fact, it can be easily checked that 

in'o", v'o') - dx = ^(-; - i)'+-V2r _ i _ ^, -i _ , 
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where T{a, x) denotes the incomplete Gamma function defined by r(a, x) = 
e~*t"'~^ dt. The previous formula can be derived by means of the change 
of variable bx^ = 2t and using [39, Formula 8.353.3]: 

T(a,x) = — 7 / dt, Rea<l, a;>0. 

' ' T{l-a) Jo t + x ' 

As wc can see, the norm of the final eigenfunction depends on I, not on 
b, and it takes real values only if I < —1, in agreement with the range of 
application for / previously derived. 



7.2. Radial Coulomb-like potentials 

Let us consider now the family of potentials 

X X'^ 

where x G (0,oo) and q ^ 0, I are real parameters. This family shares 
several characteristics with that of (83). For a start, each member can 
be regarded as being part of a pair of shape invariant partner potentials, 
respect to the transformation law Z ^ ^ + 1 [1, 9]. Similarly as before, one 
can obtain the normalized eigenfunctions (up to a modulus one factor) [36] 



'^'^ w-yr(2Z + 2 + fc)(fc + z + i)'+2^ [kTT+i)- 

(90) 

These eigenfunctions are squarc-intcgrable only in the following circum- 
stances: for I e (—3/2, —1) and q > 0, only the eigenfunction with k = 0. 
For I G (—1,00) and q < 0, the functions (90) are normalizable for all 

/c = 0, 1, 2, The normalizable solution in the range I e (—3/2,-1), 

g > 0, goes to infinity as x tends to zero, meanwhile all the others go to zero 
as a; ^ 0. Again, the reason is the existence of different self-adjoint exten- 
sions on the different ranges, being the Hamiltonians Hi^g = —-^ + Vi^q{x) 
essentially self-adjoint if I > —3/2 and q/{l -|- 1) < 0. 

The corresponding eigenvalues to the eigenfunctions (90) for the family 
(89) are [36]: 

^^" = -OfcTTTIF' ^ = M,2,... (91) 



li q — —1 and thus I > —1 we recover the spectrum given, for example, in 
[19]. Compare also with [15, p. 389]. 
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Example 7.2. Let us use now Theorem 5.2 with two potentials of the 
family (89) with different values of I. We ask that 

Vl,q{x) - Vi_r,qW = o ~ 



where r > is to be determined below, so we can choose 7(,r(a^) = 
x/\/r{2l + 1 — r). We will transform one eigenfunction Ci '^'^Ca;) which 
satisfies 



for some A; = 0, 1, 2, ...,by using some suitable solution of the equation 

A natural idea is to choose (j)v{x) as one of the eigenfunctions Cm{x) of 
Vi^q{x) for certain integer m defined by the condition 

2 



(fc + /-r + l)2 



whose simplest solution is m = fc — r. Since m and k are non-negative 
integers, we have k > r and therefore r must be a non-negative integer 
as well. As in Example 7.1, in order to avoid singularities in the final 
potential, we have to take m = and hence k = r. Then, we transform the 
eigenfunction corresponding to the integer fc > of the potential Vi-k,q{x), 
with eigenvalue —q^ /{I + 1)^, by using the ground state of the potential 
Vi^q(x), with the same energy eigenvalue. The original eigenfunction has 
to be normalizable, so it must be Z — fc > — 1 and hence g < 0, because in 
the range Z G (—3/2,-1) there are only one normalizable eigenfunction and 
A; > 0. Consequently, I > k — l>Q, and both of the initial and intermediate 
eigenfunctions are square-integrable and go to zero as a; ^ 0. If we denote 



''^(cc) = (1/Co ''(a;))(iCo ''(a;)/(ia;, the image potential reads 



V^_,,(,) - 2(^ + = V^.^x) - ^ = 1^_,,,,,,(.) . 

Correspondingly we find, after some calculations, the final eigenfunction 

l.q/ \ ( \( '^''^{x) l.qJ-k.q, ^\ I ^ ^l-k,q 1/ (l+l) f s. 



dx " ^ ' y i + i 
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In this way we recover the original potential Vi-k,q{x) but with the coupling 
constant q scaled by the factor + 1) > 0. This scaling is also reflected 
in the final eigenfunction, which moreover has k—1 instead of k, and norm 



Example 7.3. We will consider now the following modified version of 
the potentials (89): 

where again x G (0, oo) and I, q are real parameters. The normalized 
eigenfunctions are given also by (90), and as before there exist only the 
normalizable eigenfunction ioi k = if I G (—3/2, —1), q > 0, and in the 
range I G (—1, oo), q < 0, for all /c = 0, 1, 2, ... . However, the correspond- 
ing eigenvalues for the potentials (92) are now 

^^'' = (mF-(^TTTIF' ^ = «'i'2'- 

As in previous examples, we use two members of the family (92) with 
different values of I. Following Example 7.1, we think of using Vi+i^q{x) as 
the initial potential and Vi^q{x) as the intermediate one with —3/2 < I < 
— 1. The eigenfunction of the initial potential has to be square-intcgrablc 
so we must set g < 0. This means that the intermediate potential Vi^q{x) 
will have no square-intcgrablc eigenfunctions. 

One simple way to overcome this difficulty is just to change the sign of 
q in the intermediate potential, which is what we will do in this Example. 
The interesting point, however, is that it is even possible to use a non 
normalizable eigenfunction of Vi_q(x) as the intermediate one, leading to 
physically interesting results. We will see this in the next Example. From 
the analysis of these two examples it can be shown that the range I G 
(—3/2,-1) is indeed the only possibility if we restrict q to take the same 
absolute value in the initial and intermediate potentials. 

Now, assuming that q < 0,we calculate the difference 



Vl-q{x) - Vl+l,q{x) 



q' q' 2{l + l) Aq 



{l + lf (/ + 2)2 a;2 X 

with e\'^ > 0. The third and fo 
always positive for a; € (0, oo) if ^ < —1 and g < 0. An appropriate 7(0;) is 



The first two terms coincide with E^'^ > 0. The third and fourth are 
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therefore 



v^ ,%t:^;. -4gx- 2(1+1) 



The spectra (93) of two members of the family (92) with values of I differing 
by one coincide only for the ground state energy. As E^'^ — for all I, q, we 
will transform the ground state of Vi+i^q{x) by using the ground state of 
VL-q{x), both of them with zero energy. The final potential is, after some 
calculations, 

_2q (/ + l)(/ + 2) g2 



+ 



X ' x2 (Z + 2)2 

2{l + l)q{2{l + 1)(^ + 2)-^^ + {2P + 61 + 5)qx} 
^ 2(1 + If {I + 2)2(/ + 1 + 2qx)x - [21 + 2,)qH^ 
4(l + l)2(; + 2)2(2? + 3)(z3x2 



x{2{l + 1)2(/ + 2)2(1 + 1 + 2qx) - {21 + 3)g2a;2}2 
2{l + If {I + 2)2g{(2;3 + 10/2 + 10/ _ i-)qx + A{1 + if {I + 2)2} 
x{2{l + If {I + 2f{l + 1 + 29a;) - {21 + 3)92a;2}2 

where Vq~'^{x) = {l/(^Q~'^{x))d(^Q~^{x)/dx, as usual. The known eigen- 
function, with zero energy, of the previous potential is 

_ 2'+l|g|'+V2ema;'+2{(/ + !)(/ + 2) + (2/ + :\)q.r} 

{1 + 1){1 + 2)'+4 Vr(2l+4)/ ^(^'+^jgg, - MX - 2{l + 1) ■ 

Since I € (—3/2,-1) and g < 0, this function has neither zeros nor singu- 
larities in (0, 00). Moreover it is square-integrable, for the integral 



Jo 



10 

becomes after the change of variable t = 2\q\x/{l + 2), 
1 



2{l + 2)r{2l + 4) 



{A{1 + l)h{l) - 4{l + 1){21 + 3)/2(0 + {21 + 3)2/3(0} ., 
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where 

andd{l,t) = {?, + 2l)t^+%{l + 2){l + lft-8{l + lf. These integrals converge 
when I € (—3/2, —1). We have computed numerically the complete expres- 
sion and checked that it takes positive real values in the same interval. The 
result is a function strictly increasing with I, varying from approximately 
0.4 to 1. Taking into account these properties, the eigenfunction riQ^{x) 
should be the ground-state of the image potential. 



Example 7.4. As our final example we will consider the previous 
one but using a non square-integrable eigenfunction, but without zeros, 

of the intermediate potential. As sometimes happens for the standard 
intertwining technique, we will arrive to a physically meaningful image 
potential (see, for example, [4, 5]). 

Consider again the family of potentials (92). We choose Vi+i,q{x) as the 
original potential, with I e (—3/2,-1) and q < 0. The potential Vi^q{x) 
will be the intermediate one. Their associated spectra coincide just at 
zero energy, although the corresponding eigenfunction for the intermediate 
potential is not square-integrable. If we consider the difference 

T/ ^ ^ T/ ( ^ 2(/ + l) 

V^„{x) - Vi^r„{x) = - - , 

we see again that the first two terms coincide with ijj''^ > and that the 
third one is always positive for a; G (0, oo) if ? < — 1, so we can define 



Now, we transform the ground state of Vi+i^q{x) by using the formal 
mathematical eigenfunction of Vi_q{x) with zero eigenvalue, which is not 
normalizable and has no zeros. The final potential becomes now 



Vi';^ix)=Vi+r,q{x)-2[^v'^ 



lUq " dx J -fl^q (/ + 2)2 

^2g ^ (I + + 2) 2{l + l)q{2{l + !)(/ + 2)2 + (2/ + 3)qx} 



+ 



X a;2 2{l + l)^{l + 2)'^x- {2l + 3)q^x^ 

6(? + l)3(; + 2)2(2? + 3)g2 
{2{l + 1)3(1 + 2)2 - (2/ + 3)g2x2}2 ' 
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where Vq'^{x) = {1/(q'^{x)) dCQ''{x)/dx. The known eigenfunction with 
zero energy for the image potential is of the form 



dx 

^ 2'+^\q\'-''^-^c^.r'--{(l + + 2) - q.v} 

As I € (—3/2,-1) and q < 0, Vo'^ix) has no singularities for x € (0, oo) 
but has a zero at the value xq = {I + + 2)/q > 0. This function is 
square-integrable, since the integral 



/>oo 

Jo 



becomes after the change of variable t = 2\q\x/{l + 2) 

WTWWTT)^'^' + 1)^/,(0 + 4(1 + !)/,(/) + /3(0} , 

where 

roc ^-t£2l+Z+k 

m = l ^^^dt, k = i,2,3, 

and now d{l,t) = (3 + 2l)t^ — 8{l + l)"^. These integrals can be computed 
explicitly with the aid of [39, Formula 8.389.6]: 

['^^dt = ^M^^-i{e'W+(^-)V^)r(l - .,i/3^) 

Re/3>0, Re/x>0, Rev>-1. 

In our case, ^ = 1 > 0, /? = \/—S{l + l)^/{2l + 3) is real and positive for 
I G (—3/2, —1) and i' is alternatively 21 + 4, 21 + 5 and 21 + 6, all of them 
greater than —1. The final expression for Ivli'^ {^)\'^ dx is 

4{l + + 8{l + !)(/ + 2)i2{l) + 2{2l + 5)(/ + 2)13(1) 

4{l + 2){2l + 5) 

where 

ik{l) = /3(Z)2'+2+fe{e'90.fc)r(-2; - 2 - k,il3{l)) 

+e-'s('''=)r(-2; -2-k, -iP{l))} , fc = 1, 2, 3 , 
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being g{l,k) = /3(/) + (2Z+2+fc)7r/2 and /3(/) = v'-8(^ + 1)V(2' + 3). This 
function is real, positive and strictly decreasing from approximately 3 to 1 
with I e (—3/2,-1). Then, the calculated cigcnstatc should correspond to 
the first excited state of the final potential. This implies that there should 
exist a ground state eigenfunction with negative energy eigenvalue. 



8. CONCLUSIONS AND OUTLOOK 

Along this article we have studied several important facts concerning the 
problem of A-related or intertwined Hamiltonians. To this end, we have 
used simple but powerful group theoretical ideas. 

The most important results are the the following. Firstly, we have estab- 
lished the relationship between the finite difference algorithm used in [4] 
and the transformation group on the set of Riccati equations considered in 
[27] , and we have shown that the first is a particular instance of the second. 

Secondly, we have identified the group elements preserving a subset of 
Riccati equations obtained from the set of Schrodinger-like equations by 
means of the reduction procedure explained in Section 3. In this way we 
have generalized the results of the finite difference algorithm to a situation 
which seems to be new. 

Thirdly, we have approached the problem of ^-related or intertwined 
Hamiltonians in terms of the transformation group on the set of Riccati 
equations and the reduction method of Section 3, giving new insight into 
the nature of the problem. 

Finally, we have illustrated by means of some examples the use of the new 
general Theorems found at Section 5, thus generating potentials for which 
one eigenfunction and its corresponding eigenvalue are exactly known. As 
far as we know, some of these potentials have not been considered in the 
literature until now. 

These results induce new interesting questions. For instance, in the case 
of not constant 7(a;), we do not know whether it is possible to factorize 
somehow the initial, intermediate and final Hamiltonians. 

Another point concerns new applications of Theorem 5.2. On the one 
hand, it is an obvious idea to try to use potentials different from the 
oscillator-like and Coulomb-like potentials treated here. On the other hand, 
inspired by the Example 7.4, it is natural to consider explicitly other non 
square-integrable eigenfimctions, without zeros, of the intermediate po- 
tential, and see whether they lead to physically meaningful final poten- 
tials. This will involve to find the general solution of the intermediate 
Schrodinger-like equation and to restrict then the discussion to the partic- 
ular solutions without zeros in the interval of interest, even if they have 
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no physical interpretation. This means, in some sense, to adapt to our 
current method the idea introduced by Mielnik in [40], and developed later 
in [41, 42, 43, 44], among other articles. 

Finally, a finite difference formula has been used by Adler in order to dis- 
cuss the Backlund transformations of the Painleve equations [34, 35] , also 
related with what arc called dressing chains and the well-known Kortcwcg 
de Vries equation, see [45, 46] and references therein. Moreover, the Dar- 
boux transformation can be generalized by using more than one interme- 
diate cigcnfunction of the original problem [47]. Likewise, there exist gen- 
eralizations of the usual intertwining technique to spaces with dimension 
greater than one [48, 49, 50]. The natural question is whether there is some 
relationship between these topics and the transformation group on the set 
of Riccati equations discussed in this paper. 

We hope to provide some answers to these questions in future articles. 
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